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A model for J/^ suppression at a high energy heavy ion collision is presented. The main 
(and the only) reason for the suppression is J/ty inelastic scattering within hadron matter. The 
hadron matter is in the form of an evolving multi-component point-like non-interacting gas. 
The estimations of the sound velocity in the gas and the approximate pattern of cooling caused 
by the longitudinal expansion are also presented. It is shown that under such circumstances 
and with J/ty disintegration in nuclear matter added, J/fy suppression evaluated agrees well 
with NA38 and NA50 data. Additionally, some rough estimations of J /\P suppression in hadron 
matter described in the framework of the excluded volume gas model are given. 
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I. INTRODUCTION 



The existence of the so-called quark-gluon plasma (QGP) has been the one of the most intriguing questions of 
high energy physics for the almost two past decades. This novel state of matter has been predicted upon lattice 
QCD calculations (for a review see [|| and references therein) and the critical temperature T c for the ordinary 
hadron matter-QGP phase transition has been estimated in the range of 150 - 270 MeV (this corresponds to 
the broad range of the critical energy density e c ~ 0.26 — 5.5 GeV/fm 3 ). Since NA50 Collaboration estimates 
for the energy density obtained in the central rapidity region (CRR) give the value of 3.5 GeV/fm 3 for the most 
central Pb-Pb data point Q , it is argued that the region of the existence of the QGP has been reached in Pb-Pb 
collisions at CERN SPS. The main argument for the QGP creation during Pb-Pb collisions at the CERN SPS 
is the observation of the suppression of J/vf relative yield. J/*f? suppression as a signal for the QGP formation 
was originally proposed by Matsui and Satz || . The clue point of the NA50 Collaboration paper || is the figure 

(denoted as Fig. 6 there), where experimental data for Pb-Pb collision values of — (the ratio of the J/W 
to the Drell-Yan production cross-section times the branching ratio of the J /\& into a muon pair) are presented 
together with some conventional predictions. Here, "conventional" means that J/ 1 ^ suppression is due to J/^f 
absorption in ordinary hadron matter. Since all those conventional curves saturate at high transverse energy 
Et, but the experimental data fall from Et — 90 GeV much lower and this behaviour could be reproduced on 
the base of J/ty disintegration in QGP ||, it is argued that QGP had to appear during most central Pb-Pb 
collisions. In general, the immediate reservation about such reasoning is that besides those already known (see 
e.g. J6|— |l~0t| and references [12-15] in ||), there could be a huge number of different conventional models, so until 
this subject is cleared up completely, no one is legitimate to ruled out the absorption picture of J suppression. 

In the following talk, the more systematic and general description of Jl^> absorption in the framework of 
statistical analysis will be presented. The main features of the model are [[jj : 

• 1. a multi-component non-interacting hadron gas appears in the CRR instead of the QGP; 

• 2. the gas expands longitudinally and transversely; 

• 3. J /'J/ suppression is the result of inelastic scattering on constituents of the gas and on nucleons of 
colliding ions. 



"Talk based on the work done in collaboration with L.TurkoQ. 
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II. THE TIMETABLE OF EVENTS IN THE CRR 

For a given A-B collision t = is fixed at the moment of the maximal overlap of the nuclei (for more details 
see e.g. Q). As nuclei pass each other charmonium states are produced as the result of gluon fusion. After half 
of the time the nuclei need to cross each other (t ~ 0.5 fm), matter appears in the CRR. It is assumed that 
the matter thermalizes almost immediately and the moment of thermalization, to, is estimated at about 1 fm 
P,pl|. Then the matter begins its expansion and cooling and after reaching the freeze-out temperature, 2/. ., it 
ceases as a thermodynamical system. The moment when the temperature has decreased to T/ . is denoted as 
tf. .. Since the matter under consideration is the hadron gas, any phase transition does not take place during 
cooling here. 

For the description of the evolution of the matter, relativistic hydrodynamic is explored. The longitudinal 
component of the solution of hydrodynamic equations (the exact analytic solution for an (l+l)-dimensional 
case) reads (for details see e.g. puO) 



, > s T Q n° B T Q z 

s(t) = , n B (T) = -2— , v z = - (1) 

T T t 



where r = y/t 2 — z 2 is a local proper time, v z is the ^-component of fluid velocity (z is a collision axis) and sq 
and n° B are initial densities of the entropy and the baryon number respectively. For ns — and the uniform 
initial temperature distribution with the sharp edge at the border established by nuclei radii, the full solution 
of (3+l)-dimensional hydrodynamic equations is known p3[ . The evolution derived is the decomposition of the 
longitudinal expansion inside a slice bordered by the front of the rarefaction wave and the transverse expansion 
which is superimposed outside of the wave. Because small but nevertheless non-zero baryon number densities 
are considered here, the above-mentioned description of the evolution has to be treated as an assumption in the 



presented model. The rarefaction wave moves radially inward with a sound velocity c s (see Sect. IV) 



III. THE MULTI-COMPONENT HADRON GAS 

For an ideal hadron gas in thermal and chemical equilibrium, which consists of I species of particles (here, 
mesons are up to if| and baryons up to energy density e, baryon number density ub, strangeness density 
ns and entropy density s read (h = c = 1 always) 

e =_LV(2 Sl + l) / dp , (2a) 



-LY(2s 1 + 1) [ dp 

2 - 2 k Jo P exp{^} 



1 T-^ f I R 



113 O-rr 

-" i= i -"J exp 



113 2** 



6tt 2 T 2 



1 1 f°° n 2c t- 

£(2 Si + l)/ dp , (2c) 

n U J ° exp{^}+ ffl 

« roc 4 (Ei - Pi) exp j ^hp± \ 

J ° E% (exp{^}+ ff4 ) 



where Ei = (m 2 + p 2 ) 1 / 2 and m^, Bi, Si, pi, Si and gi are the mass, baryon number, strangeness, chemical 
potential, spin and a statistical factor of specie i respectively (an antiparticle is treated as a different specie) . 
And pi = Bips + Sips, where ps and ps a^e overall baryon number and strangeness chemical potentials 
respectively. 

To obtain the time dependence of temperature and baryon number and strangeness chemical potentials one 



has to solve numerically equations (2b - 2d) with s, ub and ns given as time dependent quantities. For s(r) and 
Ub(t) expressions ([j]) are taken and 715 = since the overall strangeness equals zero during all the evolution 
(see Sect. 0). 
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IV. THE SOUND VELOCITY IN THE MULTI-COMPONENT HADRON GAS 



In the hadron gas the sound velocity squared is given by the standard expression 

Since the experimental data for heavy-ion collisions suggests that the baryon number density is non-zero in the 
CRR at AGS and SPS energies [|l4]4l6|l, we calculate the above derivative for various values of n B |p7| , |l8| 



To estimate initial baryon number density rig we can use experimental results for S-S [jlij or Au-Au |l5| . [L6[ 
collisions. In the first approximation we can assume that the baryon multiplicity per unit rapidity in the CRR is 



proportional to the number of participating nucleons. For a sulphur-sulphur collision we have dN B /dy = 6 14 
and 64 participating nucleons. For the central collision of lead nuclei we can estimate the number of participating 
nucleons at 2 A = 416, so we have dN B /dy = 39. Having taken the initial volume in the CRR equal to ttR\ ■ 1 
fm, we arrive at n B = 0.25 far 3 . This is some underestimation because the S-S collisions were at a beam energy 
of 200 GeV/nucleon, but Pb-Pb at 158 GeV/nucleon. From the Au-Au data extrapolation one can estimate 
n° B = 0.65 fm -3 p5| . These values are for central collisions. So, we estimate (||) for n B = 0.25, 0.65 far 3 and 
additionally, to investigate the dependence on n B much carefully, for n B = 0.05 far 3 . The results of numerical 
evaluations of (|^) are presented in Fig.[l]. For comparison, we drew also two additional curves: for n B = and 
for a pure massive pion gas. These curves are taken from 




FIG. 1. Dependence of the sound velocity squared on temperature for various values of ub'- ns ~ 0.65 fm 
(short-dashed), ns = 0.25 fm -3 (dashed), ns = 0.05 fm -3 (solid) and ns = (long-dashed). The case of the pure pion 
gas (extra-long-dashed) is also presented. 



The "physical region" lies between more or less 100 and 200 MeV on the temperature axis. This is because 
the critical temperature for the possible QGP-hadronic matter transition is of the order of 200 MeV ||] and 
the freeze-out temperature should not be lower than 100 MeV p|. In low temperatures we can see completely 
different behaviours of cases with n B = and n B ^ 0. We think that this is caused by the fact that for n B ^ 
the gas density can not reach zero when T — ► 0, whereas for n B = it can. For the higher temperatures all 
curves excluding the pion case behave in the same way qualitatively. From T w 70 MeV they decrease to their 
minima (for n B = 0.65 far 3 at T = 219.6 MeV, for n B = 0.25 far 3 at T = 202.5 MeV, for n B = 0.05 far 3 at 
T = 183.4 MeV and for n B = at T = 177.3 MeV) and then they increase to cover each other above T w 250 
MeV. 

To compare with the more realistic picture, we also calculated the sound velocity for the excluded volume 
multi-component hadron gas model (for details of the model see [^9|). The gas consists of the same species of 



particles as in Sect. [II. The radii of all particles are assumed to be the same and equal to 0.4 fm [£(J. The 
results are presented in Fig.||. The most crucial new feature of the sound velocity is that it can be greater 
than 1 in some wide intervals of temperature. This has already been noticed [ pT| , but for the case with baryon 
number and strangeness chemical potentials equal to zero. Note that now there is strong dependence of the 
sound velocity on the baryon number density, only cases with n B = 0.0, 0.05 far 3 have c 2 s values comparable 
with the case of the point-like gas in the "physical region" (cf. Fig.[l] and Fig. ||). 
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FIG. 2. Dependence of the sound velocity squared on temperature for the excluded volume hadron gas model for 
various values of ub'- nj = 0.65 fm -3 (short-dashed), ns = 0.25 fm -3 (dashed), ns = 0.05 fm~ 3 (solid) and ns = 
(long-dashed). 



V. THE PATTERN OF COOLING AND ITS CONNECTION WITH THE SOUND VELOCITY 



In Sect. HI we have explained how to obtain the time dependence of the temperature of the longitudinally 
expanding hadron gas. This dependence proved to be very well approximated by the expression |I^ , |l"8f| 



T(t) = T -t 



(4) 



The above approximation is valid in the temperature ranges [T/. ., To], where T/. . > 100 MeV, To < To.max and 
T , ma x « 230 MeV. We started from T equal to 227.2 MeV (for n% = 0.65 fm~ 3 ), 229.3 MeV (for n% = 0.25 
fm -3 ) and 229.7 MeV (for n B = 0.05 fm~ 3 ). These values correspond to eo = 5.0 GeV/fm 3 (the initial energy 
density in the CRR has been estimated by NA50 at eo = 3.5 GeV/fm 3 ). Then we took several decreasing 
values of To < T^, max . For every To chosen we repeat the procedure of obtaining the approximation (^), i.e. 
the power a. The values of a as a function of To are depicted in Fig. || for the case of tib = 0.25 fm~ 3 , together 
with the corresponding sound velocity curve. We can see that in the above-mentioned interval of temperature 
the power a has the meaning of the sound velocity squared within quite well accuracy. We would like to stress 
that we arrived at the same conclusion also for cases with ub — 0.65 fm -3 and ns — 0.05 fm -3 . 



1 00 200 

temperature [MeV] 



FIG. 3. The power a (crosses) from the approximation (Ul) and the sound velocity curve for tib — 0.25 fm 
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We can formulate the following conclusion: in the "physical region" of temperature and for realistic baryon 
number densities, the longitudinal expansion given by ([!]) results in the cooling of the hadron gas described by 
(g) with a = c^(T ), namely: 

T(r) = T .(^) (5) 

2 

where Tq belongs to the "physical region". Note that T(r) = T ■ (^) Cs is the exact expression for a baryonless 
gas with the sound velocity constant (for details see |^2|,^3|). It should be stressed that c 2 s in (|5|) is constant 
only for the particular cooling from a given To. For different To the value of c 2 s in (||) will change appropriately, 
that is to the value of c 2 s in the hadron gas with the temperature equal to this new T . The approximation (j^) 
will be used to simplify evaluations of J/\I / survival factors in next sections. 

It should be added that for the excluded volume hadron gas model the above-mentioned conclusion is no 
longer valid. We did the appropriate simulations, but approximations of T(t) by a power function were much 
worse and values of the power obtained were a factor 2 lower than the square of the velocity of sound. 



VI. J/vp ABSORPTION IN THE EXPANDING HADRONIC GAS 

As it has been already mentioned in Sect. |^ charmonium states are produced in the beginning of the collision, 
when nuclei overlap. For simplicity, it is assumed that the production of cc states takes place at t = 0. To 
describe J/^> absorption quantitatively, the idea of || is generalized to the case of the multi-component massive 
gas, here. Since in the CRR longitudinal momenta of particles are much lower than transverse ones (in the 
c.m.s. frame of nuclei), J/^ longitudinal momentum is put at zero. Additionally, only the plane z = is under 
consideration, here. For the simplicity of the model, it is assumed that all charmonium states are completely 
formed and can be absorbed by the constituents of a surrounding medium from the moment of creation. The 
absorption is the result of a cc state inelastic scattering on the constituents of the hadron gas through interactions 
of the type 

cc+h — >D + D + X, (6) 

where h denotes a hadron, D is a charm meson and X means a particle which is necessary to conserve the 
charge, baryon number or strangeness. 

Since the most crucial for the problem of QGP existence are Pb-Pb collisions (see remarks in Sect. the 
further considerations are done for this case. So, for the Pb-Pb collision at impact parameter b, the situation 
in the plane z — is presented in Fig.^J, where S e /f means the area of the overlap of the colliding nuclei. 



FIG. 4. View of a Pb-Pb collision at impact parameter b in the transverse plane (z = 0). The region where the nuclei 
overlap is hatched and denoted by S e //. 

Additionally, it is assumed that the hadron gas, which appears in the space between the nuclei after they 
have crossed each other, also has the shape of S e /f at to in the plane z — 0. Then, the transverse expansion 
starts in the form of the rarefaction wave moving inward S e ft at to. 

From the considerations based on the relativistic kinetic equation (for details see p],[l ) , the survival fraction 
of J/vff in the hadron gas as a function of the initial energy density eo in the CRR is obtained: 
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A^. 9 .( £ o) = / dp T g(pr, eo) ■ exp 



to 



(2tt 



sfi{q,t)aiV re i ti ^^ 



(7) 



where the sum is over all taken species of scatters (hadrons), p v — (E,px) and q\ — (E' i7 q) are four momenta 
of J /if? and hadron specie i respectively, Oi states for the absorption cross-section of J/^ — hi scattering, u r ei,i 
is the relative velocity of hi hadron with respect to J /if? and M and rrii denote J /if? and hi masses respectively 
(M = 3097 MeV). The function g(pT,eo) is the J /if? initial momentum distribution. It has a gaussian form 
and reflects gluon multiple elastic scatterings on nucleons before their fusion into a J /if? in the first stage of 
the collision [p2|-p5|. The upper limit of the integration over time in (Q), namely t final is the minimal value 
of (t esc ) and i/. . • The quantity (t esc ) is the average time of the escape of J/^f's from the hadron medium for 
given values of b and J /if? velocity v = pr/E. Note that the average is taken with the weight 



pj/*{r) 



T A (f)T B (r-b) 
T AB (b) 



(8) 



where TabQ>) = J cPs ' Ta(s)T b (s — b), Ta(s) = J dzpA{s, z) and pa(s,z) is the nuclear matter density distri- 
bution. For the last quantity, the Woods-Saxon nuclear matter density distribution with parameters from (20) 
is taken. In the integration over hadron momentum in (Q) the threshold for the reaction (g) is included, i.e. <Ji 
equals zero for (jp v + q 1 /) 2 < (2m B + nix) 2 and is constant elsewhere (m B is a charm meson mass, mo = 1867 
MeV). Also the usual Bose-Einstein or Fermi-Dirac distribution for hadron specie i is used in (Q) 



fi(q,t) = fi(q,t) 



2s, + 1 



exp 



(9) 



For simplicity, we use (||) as the approximation to T(t) in (||) and p B (t) and ps(t) are solutions of only two 
equations @ - 1§) with T given by (g), n B (t) by ([[]) and n s (t) = 0. 

As far as o~i is concerned, there are no data for every particular J /if 1 — hi scattering. Therefore, only two types 
of the cross-section, the first, Ob, for J /if? — baryon scattering and the second, er m , for J /if? — meson scattering 



are assumed, here. In the model presented <j\> is put at a 



Nucleon absorption cross-section. 
The J /if? — meson cross-section 



And values of oj/^n in the range of 3-5 mb are taken from p-A data 
<r m is assumed to be 2/3 of o~b, which is due to the quark counting. 

As it has been already suggested |27]] also J /if? scattering in nuclear matter should be included in any J /if? 
absorption model. This could be done with the introduction of J/if? survival factor in nuclear matter (29^33] 



A/" n . OT .(eo) = exp {-o-j/^nPqL} 



(10) 



where po is the nuclear matter density and L the mean path length of the J/if? through the colliding nuclei. 
The length L is expressed by the following formula 132] : 



PoL(b) 



1 



2Tab 



d 2 sT A {S)T B {s~b) 



A-l 
A 



B — 1 
B 



T B (s-b) 



(11) 



Since J/if? absorptions: in nuclear matter and in the hadron gas, are separate in time, J/if? survival factor 
for a heavy-ion collision with the initial energy density e , could be defined as 

Af(e ) = Af n .m.{e ) ■ Af h .g.(eo) ■ (12) 

Note that since right sides of (0) and ([[o|) include parts which depend on impact parameter b and the left sides 
are functions of eo only, the expression converting the first quantity to the second (or reverse) should be given. 
This is done with the use of the dependence of eo on the transverse energy Et extracted from NA50 data [| 
(for details see 

To make the model as much realistic as possible, one should keep in mind that only about 60% of J/if?'s 
measured are directly produced during collision. The rest is the result of \ (~ 30%) and tp' (~ 10%) decay [ |33[ . 
Therefore the realistic J/if? survival factor could be expressed as 

Af(eo) = 0.6AA w (e ) + 0.3A/" x (e ) + 0.LV^(e o ) , (13) 

where A/j/v>( £ o)j N x (eo) and A^'(eo) are given also by formulae (Q), ( [l0| ) and ( |l2] ) but with g[pT,eo), &j/ipN 
and M changed appropriately (for details see [0). 

To complete, also values of cross-sections for \~ baryon and tp' — baryon scatterings are needed. For simplicity, 
it is assumed that both these cross-sections are equal to J /if? — baryon one. Since J /if? is smaller than \ or -ip' , 
this means that J /if? suppression evaluated according to (O) is underestimated here. 
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VII. RESULTS 



Now we can complete calculations of formula (0) for values of n° B given in Sect. [V. The last parameter of 
the model is the freeze-out temperature. Two values T/. . = 100, 140 MeV are taken here and they agree well 
with estimates based on hadron yields |lq] . 

Firstly, formula ( p"3| ) with the use of (]?]) only instead of ( |l2| ) is calculated for the case of (l+l)-dimensional 
expansion and results are presented in Figs.|^-p[ 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

£ [GeV/fm 3 ] 
FIG. 5. J/^ suppression in the longitudinally 
expanding hadron gas for n% — 0.25 fm -3 and 
Tf. . = 140 MeV: (a) o b = 3 mb, a m = 2 mb; (b) 
a b = 4 mb, a m = 2.66 mb; (c) a b = 5 mb, a m = 3.33 
mb; (d) a b = 6 mb, <r m = 4 mb . The black squares 
correspond to the NA38 S-U data Q, the black tri- 
angles correspond to the 1996 NA50 Pb-Pb data, the 
white squares to the 1996 analysis with minimum bias 
and the black points to the 1998 analysis with mini- 
mum bias 




1.5 2.0 2.5 3.0 3.5 
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0.0 0.5 1.0 

£ iGeV/fm 
FIG. 6. Same as Fig.H (except case (d), which is 
not presented here) but for Tf.„. = 100 MeV. 



For comparison, also the experimental data are shown in Figs. 



The experimental survival factor is defined 



as 



Nt ?.xp — 



' DY 



"TTTO' 



(14) 



wnere "abVpp) — i s t ne ratio of the J/^P to the Drell-Yan production cross-section in A-B(p-p) interactions 

DY 

times the branching ratio of the J/vE' into a muon pair. The values of the ratio for p-p, S-U and Pb-Pb are 
taken from |3 34 3^] . The results of numerical evaluations of ( |T5| ) agree well with the data for greater baryonic 
cross-section ct, and (or) for the lower freeze-out temperature T/. ., except the last point measured for the 
highest Et bin. 

Now the full (3+l)-dimensional hydrodynamic evolution (in the form described in Sect. and J/'F absorption 
in nuclear matter will be taken into account and results are depicted in Figs.fj||| To draw also S-U data together 
with Pb-Pb ones, instead of multiplying Afh.g. 
hadron gas survival factor" is defined as 



by A/'n.m., Nexp is divided by Af n .m., i-e. "the experimental J/^ 



Nexp = exp {{Jj/^nPqL} ■ Nexp 



(15) 



and values of this factor are drawn in Figs.^-|| as the experimental data. The quantity r is a parameter from 
the expression for a nucleus radius Ra — fo • -As. The value of Ra is used to fix the initial position of the 
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rarefaction wave in the evaluation of (t esc ) in (M) 
differ substantially from that of n° B = 0.25 fm 



It has turned out that the case of n B 



so curves for n B 



0.05 far 



0.05 fm 3 does not 
are not depicted in Figs.^-^. 




0.0 0.5 




FIG 



£ [GeV/frri ] 



W 

8. Same as Fig. 171 but for at 
3.33 mb. 



= 5 mb and 



1.0 1.5 2.0 2.5 3.0 3.5 

£ [GeV/fm 3 ] 
FIG. 7. J/\& suppression in the longitudinally 
and transversely expanding hadron gas for the 
Woods-Saxon nuclear matter density distribution and 
a b = 4 mb, a m = 2.66 mb and T f . . = 140 MeV. The 
curves correspond to n° B = 0.25 fm~ 3 , c B — 0.45, 
r = 1.2 fm (solid), n% = 0.65 fm~ 3 , c 3 = 0.46, 
r = 1.2 fm (dashed) and n% = 0.25 fm -3 , c s = 0.45, 
ro = 1.12 fm (short-dashed). The black squares rep- 
resent the NA38 S-U data |34j , the black triangles rep- 
resent the 1996 NA50 Pb-Pb data, the white squares 
the 1996 analysis with minimum bias and the black 
points the 1998 analysis with minimum bias ||, but 
the data are "cleaned out" from the contribution of 
J /\P scattering in the nuclear matter in accordance 
with @. 

Having compared Figs. with Figs.|[]^, one can see that adding the transverse expansion changes the hnal 
(theoretical) pattern of J/\f f suppression qualitatively. Now the curves for the case including the transverse 
expansion are not convex, in opposite to the case with the longitudinal expansion only, where the curves are. 
As far as the pattern of suppression is concerned, theoretical curves do not fall steep enough at high eo to 
cover the data area. But for some choice of parameters, namely for ab somewhere between 4 and 5 mb and 
for ro = 1.12 fm, a quite satisfactory curve could have been obtained. Precisely, again only the highest Er 
bin point falls down of the range of theoretical estimates completely. But the error bar of this point is very 
wide. And also the contradiction in positions of the last three points of the 1998 data can be seen. This means 
that the high Et region should be measured once more with the better accuracy to state definitely whether the 
abrupt fall of the experimental survival factor takes place or not there. 

To support the conclusion, main results from Figs.0-^ are presented in Fig.|| again. The original data j|] for 

B a PbPb B a VV 

— PfrPif and J/\I> survival factors given by (U) multiplied by "V^ and now as functions of Et are depicted 
there (for details see fm). 
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FIG. 9. J/ty survival factor times pp^ in the longitudinally and transversely expanding hadron gas for the 

a r>Y 

Woods-Saxon nuclear matter density distribution and n° B = 0.25 fm -3 , Tf, , = 140 MeV, c s — 0.45 and ro = 1.2 fm. The 
curves correspond to Ob — 4 mb (solid) and at = 5 mb (dashed). The black triangles represent the 1996 NA50 Pb-Pb 
data, the white squares the 1996 analysis with minimum bias and the black points the 1998 analysis with minimum bias 

1- 



As it has been already mentioned, the main disagreement with the data appears in the last experimental point 
of the 1998 analysis. 

The charmonium — baryon inelastic cross-section is the most crucial parameter in the model presented. So 
to be sure what is the exact result of J/\l/ absorption, one should know how this cross-section behaves in the 
hot hadron environment. The newest estimations of tt + J/^>, p + J/^ and J/^! + N cross-sections at high 
invariant collision energies |36 37]] give values of ah and a m of the same order as assumed here. Also it should 
be stressed that the charmonium — hadron inelastic cross-sections are considered as constant quantities in 
the model presented. But, as the results of just mentioned papers suggest, they are not constant at all. The 
cross-sections are growing functions of the invariant collision energy \fs. Therefore, one could think naively 
that the increase of eo (or in other words Et) causes the increase of the invariant collision energy ^fs on the 
average and further the increase of the charmonium — hadron inelastic cross-sections. In this way, the line 
describing J/vP survival factor could be placed close to the solid curve of Fig. ^| for low e {Et), and then, as the 
charmonium — hadron inelastic cross-sections increase, this line would go closer to the dashed curve of Fig.|| 
for high £q {Et)- So, the experimental pattern of J/\P suppression could be recovered. 



VIII. THE SUPPRESSION IN THE HADRON MATTER DESCRIBED BY THE EXCLUDED 

VOLUME GAS MODEL 

As a last remark an extension to the point-like gas model will be discused. The rough estimation of J /if? 
suppression in the framework of the excluded volume gas model Jl9[ will be presented now. The estimation done 
suggests that J /if? suppression is much lower in the gas of particles with hard core radii than in the point-like 
one. There are two main reasons for such behaviour: 

• 1. the freeze-out time resulting from the longitudinal expansion only, is much shorter for the excluded 
volume gas than for the point-like one, namely for n? B = 0.25 fm -3 , T/. D . = 140 MeV and eo = 3.5 
GeV/fm 3 one obtains tj x v ' — 9.4 fm and t^'~ 1 ~ 15.1 fm {"ex.v." means "excluded volume" and "p. — /." 
"point-like"); 

• 2. the initial particle density is much lower for the excluded volume gas than for the point-like one, that 
is nff = 0.8 fm~ 3 and n p ^ L = 2.6 fm" 3 . 

To estimate how these two facts influence J/vE' survival factor the very rough approximation to (j/J) is explored. 
If the transverse expansion, particle masses and thresholds are neglected, then (0) takes the following very simple 
form: 
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Mh.g. = cxp|-(Tn toln^^| , (16) 

where a means one common cross-section for all particles. So, for values of , t e j?^' and n^'~ L , riQ X ' v ' given 
above and for a — o~ m — 3.33 nib one obtains: 

tq-- 1 - ~ 0.1 , Nh X g v ' * 0.6 . (17) 

For comparison Af exp — 0.1 — 0.2 there. Therefore, J/^l survival factor evaluated for the excluded volume 
gas model is far above experimentally allowed values. This means that J/\f' suppression in a hadron gas is very 
much model dependent. 



IX. CONCLUSIONS 



To sum up and to stress how vague the idea of J/ty suppression as a signal of QGP is, the author would like 
to call reader's attention to Fig. 5 of Q. The solid and dotted curves in that figure are exactly the same as 
curves just presented in Fig.||. But the distortion in the QGP is the main reason for J/\l/ suppression in j38j! 
It is also stated there, that results shown "provide evidence for the production of the quark-gluon plasma in 
central high-energy Pb-Pb collisions" . So, one could say as well that the NA50 Pb-Pb data provides evidence 
for the production of the thcr malize d and confined hadron gas in the central rapidity region of a Pb-Pb collision. 
But estimations done in Sect. VII] have just shown that it could not be true. 

Therefore, the final conclusion is the following: J/vff suppression is not a good signal for the quark-gluon 
plasma appearance at a central heavy-ion collision. 
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